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ON THE NEUMANN PROBLEM FOR MONGE-AMPERE TYPE EQUATIONS 


FEIDA JIANG, NEIL S. TRUDINGER, AND NI XIANG 


Abstract. In this paper, we study the global regularity for regular Monge-Ampere type equations 
associated with a semilinear Neumann boundary conditions. By establishing a priori estimates for 
second order derivatives, the classical solvability of the Neumann boundary value problem is proved 
under natnral conditions. The techniques build upon the delicate and intricate treatment of the standard 
Monge-Ampere case by Lions, Trudinger and Urbas in 1986 and the recent barrier constructions and 
second derivative bounds by Jiang, Trudinger and Yang for the Dirichlet problem. We also consider 
more general oblique boundary value problems in the strictly regular case. 


1. Introduction 

In this paper, we consider the following semilinear Neumann boundary value problem for the Monge- 
Ampere type equation 

(1.1) det[D‘^u — A{x,u, Du)] = B{x,u, Du), in Q, 

(1.2) DyU = ip{x,u), on dVL, 

where 11 is a bounded domain in n dimensional Euclidean space M” with smooth boundary, Du and 
D^u denote the gradient vector and the Hessian matrix of the second order derivatives of the function 
tt : n —>■ M respectively, A is a given n x n symmetric matrix function defined on 11 x M x M"', R is a 
positive scalar valued function on 11 x M x R”, (/? is a scalar valued function defined on dQ x R and u is 
the unit inner normal vector field on dQ. As usual, we use x, z, p, r to denote points in II, R, R”, 
respectively. A solution u G (7^(11) of equation (jl.ll) is elliptic when the augmented Hessian matrix 
Mu = D^u — A{x, u, Du) is positive definite, that is Mu > 0, which implies R > 0. Also, a function u 
satisfying Mu > 0 is called an elliptic function of the equation dni). Since the matrix A determines the 
augmented Hessian matrix Mu, we also call an elliptic solution (or function) an A-admissible solution 
(or function) or, by analogy with the case A = 0, an A-convex solution (or function). 

We shall establish an existence theorem together with a priori estimates for elliptic solutions of the 
Neumann boundary value problem (fTTD-ra in this paper, which extend the special case where A 
is independent of p in HZ]. For this purpose, we need appropriate assumptions on A, R, p and Q. 
Assume that the matrix A is twice differentiable with respect to p and A, R and (p are differentiable 
with respect to z. Following |24] . we call the matrix A regular in II if A is co-dimension one convex 
with respect to p, in the sense that 

(1-3) Aij^u{x,z,p)iiijr]km > 
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for all {x,z,p) G 0 X M X R”", ^ R”, ^ -L r/, where = Dp^p^Aij. If the inequality (11.311 is strict, 

then the matrix A is called strictly regular. We also define the matrix A to be non-decreasing, (strictly 
increasing), with respect to z, if 

(1.4) DzAij{x,z,p)(iCj > 0, (> 0), 

for all {x,z,p) G n X R X R”, ^ G R"". The inhomogeneous term B and boundary function p are also 
non-decreasing, (strictly increasing), with respect to z, if 

(1.5) Bz{x,z,p) >(),{> Q), 

for all (x,z,p) G n X R X R” and 

( 1 . 6 ) ipz{x,z) >0,{> 0), 

for all {x, z) G <911 x R. Note that if we write the boundary value problem dnD-do]) in the general 
form 

(1.7) B[u] := F{x,u, Du, D'^u) = 0, in 11, 

(1.8) Q[u] := G{x,u, Du) = 0, on 911, 

where F and G are defined by 

(1.9) F{x, z,p, r) = det[r — A{x, z,p)\ — B{x, z,p), 

(1.10) G{x, z,p) = V ■ p — ip{x, z), 

then A,B and p non-decreasing, (strictly increasing), in z, correspond to the standard monotonicity 
conditions, Fz < 0, < 0, {Fz < 0, Gx < 0) for symmetric matrices r satisfying r > A{x, z,p), that is 

for points {x,z,p,r) G H x R x R"" x R”'^"', where F is elliptic. 

As with m, we also need the domain H to satisfy an appropriate uniform convexity condition. 
Adapting [23], we define the domain H to be uniformly A-convex, (A-convex), with respect to the 
boundary function p and an interval valued function X on 911 if H G G^ and 

( 1 . 11 ) {DiVj{x) - Dp^Aij{x, z,p)vk)TiTj < 0 , (< 0 ), 

for all (x, z,p) G 9H x R x R”, satisfying p ■ i>{x) > p{x, z), z G X{x) and vectors r = t(x) tangent to 
9H. For a given function uq on 911, we define H to be uniformly A-convex, (A-convex), with respect 
to if and uq if (|l.lip holds for all p ■ v{x) > (p{x, uo{x)), that is X = {rto}- 

From the regularity of A (II.3p . we can equivalently replace the boundary inequality p ■ u > p{x, z) 
by the boundary equality p ■ u = p{x,z), in the above definitions, as Dp^Aij{x, z,p)TiTj is then non¬ 
decreasing with respect to p^. This leads us to a further definition which is independent of the boundary 
condition Namely H is uniformly A-convex with respect to u G G^(ll) if 

(1.12) {DiVj - Dp^Aij{-,u,Du)vk)TiTj <-5 q, on 9H, 

for all vectors r = t{x) tangent to 9H. Accordingly if A is regular, H is uniformly A-convex with 
respect to p and u and u satisfies (|1.2p . it follows that H is uniformly A-convex with respect to u. 

In order to use the regularity of A in its most general form, we will need to assume the existence of 
a supersolution u to o satisfying 

(1.13) det[D'^u — A{x,u, Du)] < B{x,u, Du), in Q, 
together with the same boundary condition, 

(1.14) DjyU = ip{x,u), on 9n. 

We then have the following global second derivative estimate. 
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Theorem 1.1. Let u G (7^(0) n (7^(0) be an elliptic solution of the Neumann problem (in])-([L2i) in 
a domain Q C M”, which is uniformly A-convex with respect to u, where A G (7^(Q x M x M") 
is regular and non-decreasing, i3 > 0, G (7^(0 x M x M”) is non-decreasing and ip G x M) is 

non-decreasing. Suppose there exists an elliptic supersolution u G satisfying (I1.13h - (|1.14l) . Then 

we have the estimate 

(1.15) sup < (7, 

n 

where C is a constant depending on n, A, p,6o, and 

Theorem o extends Theorem 3.3 in m except for the supersolution hypothesis as a supersolution 
is constructed in m in the course of the proof. We also point out that, as in the restriction to 
the Neumann condition is critical for our proof and moreover as shown by the Pogorelov example, (see 
[3T] , [33] )) one cannot generally expect second derivative estimates and classical solutions of (in])-([L2|) 
for ^ = 0, when the geometric normal n is replaced by an oblique vector /3 satisfying ft ■ n > 0, that is 
in (flTOD . 

(1.16) G{x,z,p) = P ■p-p{x,z), 

no matter how smooth /?, p, B and dLl are. However if the matrix function A is strictly regular on 
0 so that we have a positive lower bound in (11.311 when z and p are bounded, then the proof is much 
simpler and also embraces oblique boundary conditions. Moreover in this case the monotonicity and 
supersolution hypotheses in Theorem 11.11 can be dispensed with. Typically second derivative behaviour 
for equation dni) in the strictly regular case is closer to that for uniformly elliptic equations while the 
challenge in the general case is to carry over the more intricate Monge-Ampere case, A = 0. Following 
we can also relax the supersolution hypothesis for uniformly convex domains in the special case 
when DpxA = 0, that is 

(1.17) A{x,z,p) = Ao{x,z) + Ai{z,p), 
where Aq G (7^(0 x M) and Ai G C^(]R x M") is regular. 

From Theorem 11.11 we obtain classical existence theorems for (in])-([L2|) under further hypotheses 
ensuring estimates for solutions and their gradients. For solution estimates, by virtue of the comparison 
principle we can simply assume the existence of bounded subsolutions and supersolutions. 

However more specific conditions for solution bounds will be treated in Section [3] of this paper, 
including an extension of the Bakel’man condition in Theorem 2.1 of HZ). For the gradient estimate 
we adopt the same structure condition used for the Dirichlet problem in |9|, namely 

(1.18) A{x,z,p) >-p.o{l + \pf)I, 

for all X G n, \z\ < Mq, p G M”" and some positive constant po depending on the constant Mq. Condition 
(ll.lSp provides a simple gradient bound for A-convex functions u in terms of a lower bound for Dj^u 
on the boundary. Combining the second derivative bounds with the lower order bounds and the global 
second derivative Holder estimates as in [laiiiiiiziEsi, we establish the following existence result by 
the method of continuity. 

Theorem 1.2. Suppose that A,B,p,u and Q satisfy the hypotheses of Theorem M.R with either A, B 

or p being strictly increasing. Assume also condition (I1.18P and that there exists an elliptic subsolution 

u G (7^(12) n (7^(0) of equation (|l.ll) . with DpU > p{-,u) on dLl and that 12 is uniformly A-convex 

with respect to p and Z = \u,u], in the sense of (II.lip . Then the Neumann boundary value problem 

(HID-dlll) has a unique elliptic solution u G (7^’" (12) for any a < 1. 
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The uniqueness of solutions follows from the comparison principle for elliptic solutions of general 
oblique boundary value problems, dEZD-dTSI); see Lemma EH The regularity for the solution u in 
Theorem Ol can be improved by the linear elliptic theory [5] if the data are sufficiently smooth. For 
example, if A, B, ip and <90 are C°°, then the solution u € From the monotonicity of (p, it is 

also enough to assume (II.lip only holds for p- v > ip{-,u) and u < z < u. Moreover if A is independent 
of z, there is no need for the last inequality. Also taking account of our remarks after the statement of 
Theorem 11.11 we only need to assume the supersolution u satisfies (jl.lSp at points where it is elliptic 
and the boundary inequality Di,u < ip{-,u), instead of (11.141) . if either A satisfies (I1.17P with Q also 
uniformly convex or A is strictly regular in Q. 

The regular condition was originally introduced in [22] in its strict form for interior regularity of 
potential functions in optimal transportation with the weak form (|1.3p subsequently introduced in 
|30j for global regularity; (see also [M])- It was subsequently shown to be sharp for regularity of 
potential functions in [21]. Optimal transportation equations are special cases of prescribed Jacobian 
equations, which have the general form, 

(1.19) I det DY (•, u, Du)\ = u, Du), 

where T is a mapping from x M x M”" into M"", ip is a non-negative scalar valued function on 
0 X M X M"’. Assuming detYp ^ 0, we see that for elliptic solutions, equation (I1.19p can be written in 
the form dnn with 

(1.20) A = -Y-\Y^ + Y,®p), S = (detyp)-V. 

The natural boundary value problem for the prescribed Jacobian equation is the second boundary value 
problem to prescribe the image, 

(1.21) Tu{n) := Y{-, u, Du){n) = Q*, 

where is another given domain in M"". The global regularity of the second boundary value problem 
(jl.l9p - (ll.2ip has already been studied in [H ESI ESI ESI EZ] for different forms of the mapping Y. As 
shown in |30| in the optimal transportation case and in |25] in the general case, condition (|1.21ll implies 
an oblique nonlinear boundary condition for elliptic functions u, that is m holds for a function 
G eC\dnxRx M*") with 

(1.22) Gp{-,u, Du) ■ iz > 0, on dQ. 

The crucial estimate in these papers is the control on the obliqueness, that is an estimate of the form, 
Gp ■ V > 5 for a positive constant 5 and this is done in |30] in the optimal transportation case, and 
extended to the general case in |25] . under appropriate uniform convexity conditions on the domain 
and target, with the latter equivalent to the uniform concavity of the function G with respect to the 
p variables. Because we are defining obliqueness with respect to the inner normal, in agreement with 
m, our function G is the negative of that in [321 Ea EQj. Once the obliqueness is estimated, the 
boundary second derivative bounds follow in [sniEniiii from the same uniform convexity conditions, 
together with the regular condition (|1.3p . similarly to the Monge-Ampere case in [33] • Note that the 
uniform concavity of G excludes the Neumann condition treated here and moreover the derivation 
of the boundary G^ estimate is much simpler, being somewhat analogous to using the strict regular 
condition. We also point out a recent paper [2] considering optimal transportation on a hemisphere 
where the obliqueness is estimated without using any uniform convexity of domains, which still gives 
the boundary estimate in the two dimensional case. Prescribed Jacobian equations also arise in 
geometric optics where solutions correspond to reflectors or refractors transmitting light rays from a 
source to a target with prescribed intensities; (see for example [351 la Eg ESllZl fT8] and references 
therein). 

On the geometric side, the Neumann boundary value problem in the more general context of aug¬ 
mented Hessian equations on manifolds arises in the study of the higher order Yamabe problem in 
conformal geometry; (see [a El ni la m)- To explain this we let (A4,g) be a smooth compact 
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Riemannian manifold of dimension n > 3 with nonempty smooth boundary dM, let Ag denote the 
Schouten tensor of the metric g and let X{Ag) = (Ai(^g), • • • , Xn{Ag)) denote the eigenvalues of Ag. 
Let r C M” be an open convex symmetric cone with vertex at the origin and / be a smooth symmetric 
function in L. The fully nonlinear Yamabe problem on manifolds with boundary is to find a metric g 
in the conformal class of the metric g with a prescribed function of eigenvalues of the Schouten tensor 
and prescribed mean curvature. For example, for a given constant c E M, we are interested in finding 
a metric g conformal to g such that 


(1.23) 


F{Ag) := f{\{Ag)) = 1, 

hg = C, 


X{Ag) E r on Ad, 
on clAd, 


where hg denotes the mean curvature of dM with respect to the inner normal. Writing g = e~‘^'^g for 
some smooth function u on Ad, by the transformation laws for the Schouten tensor and mean curvature, 
the problem (11.231) is equivalent to the following semilinear Neumann boundary value problem 


(1.24) 


/(A,(t/)) = e 

du 


—2u 


^ 9 '> 


Xg{U) E r on Ad, 
on clAd, 


with 


U = + du® du 


\\^u?g + A 




where Xg{U) denotes the eigenvalues of U with respect to g, v is the unit inner normal vector field 
to clAd, and V denotes the Levi-Civita connection with respect to g. If we choose / = det and F = 
Fn := {A = (Ai, • • • , An) E M"" : ^ A^ > 0}, then we have an example (|1.24p of a semilinear Neumann 
boundary value problem (ini-dLSl) for a Monge-Ampere type equation. In conclusion a prescribed 
mean curvature fully nonlinear Yamabe problem (|1.23p is equivalent to a semilinear Neumann problem 
(I1.24P for an augmented Hessian equation. The corresponding matrix functions in these cases will 
be strictly regular when expressed in terms of local coordinates so that in the Monge-Ampere case 
strong local estimates are available, with second order estimates being considerably simpler than the 
general regular case we treat here. In fact, the particular Neumann boundary value problem (11.2411 
with / = det has already been studied in m- In the special case of Euclidean space M”, the matrix A 
is given by 


(1.25) 




in which case our A-convexity condition (jl.lll) reduces to simply + </? > 0, (> 0), where ki denotes 
the minimum curvature of dfl. The overall organisation of this paper follows that of the Dirichlet 
problem case [9], where again the main issue was to deal with the general case of regular A. Also 
here the strictly regular case is considerably simpler in the case of smooth data but in the optimal 
transportation case with only Holder continuous densities local and global second derivative estimates 
were obtained in 1201 E], in agreement with the uniformly elliptic case. In Section [2] we prove Theorem 
11.11 which constitutes the heart of the paper. In Section [3] we provide the gradient estimate to complete 
the proof of Theorem 11.21 along with alternative solution bounds for more general oblique boundary 
value problems. In the optimal transportation case we also prove a Bakel’man type estimate for 
solutions which extends the Monge-Ampere case in HZ]. In Section 0] we switch to the strictly regular 
case and prove first and second derivative bounds for general oblique boundary value problems (II.8j) . 
where G is concave with respect to the p variables, which extend the semilinear conditions (|1.16l) . For 
this purpose we extend our definition of A-convexity so that a domain Q is uniformly A-convex, 
(A-convex), with respect to G and an interval Z if (|l.lll) holds for all (x, z,p) E dfi x M x satisfying 
G{x,z,p) > 0, z G Z and vectors r tangent to dXl. When G is independent of z, this corresponds to 
the c-convexity conditions from optimal transportation [241130] and more generally to the Y-convexity 
conditions for prescribed Jacobian equations in [25]. Finally we remark that a general theory of oblique 
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boundary value problems for augmented Hessian equations, which embraces our results in Section 4, is 
presented in [5]. 


2. Second derivative estimates 

In this section, we shall derive the second order derivative estimates and complete the proof of 
Theorem 11.11 by taking full advantage of the assumed supersolution u. Note that we only need to 
get an upper bound for the second derivatives, since the lower bound can be derived from the ellipticity 
condition D^u — H > 0. 

For the arguments below, we assume the functions y?, u can be smoothly extended to 11 x M and 
Q respectively. We also assume that near the boundary, is extended to be constant in the normal 
directions. From the equation we have 

(2.1) F[u] := log det[D‘^u — A{-,u, Du)] = B{-, u, Du), 

where B = logB. We have where {wij} = {uij — Aij} denotes the 

augmented Hessian matrix, and {w^^} denotes the inverse of the matrix {wij}. We now introduce the 
following linearized operators of F and (|2.1D . 

(2.2) L = {Dij — Dp^Aij{-,u,Du)Di), Cv = Lv — Dp^B{-,u,Du)Di. 

For convenience in later discussion, we denote D^^u = DijU^irjj, = Wij^irjj = DijU^itjj — Aij^irjj for 
any vectors ^ and r]. As usual, C denotes a constant depending on the known data and may change 
from line to line in the context. 

Before we start to deal with the second derivative estimates, we recall a fundamental lemma in 
which is also crucial to construct the global barrier function using the supersolution in our situation. 
We shall omit its proof, which is similar to those in mu- 

Lemma 2.1. Let u G C‘^{Cl) be an elliptic solution of (II.ip . u € C‘^{Q) be an elliptic function of 
equation dni) in Q, with u > u in Ll, where A is regular and non-decreasing. Then 

(2.3) C ~ 

i 

holds in Q for sufficiently large positive eonstant K and uniform positive constants ei,C depending on 
A, B, Ll, |u|i;n and u. 

We assume that the domain Q is uniformly A-convex, with respect to ip and u, and first consider 
the second derivative estimates on the boundary dLl in nontangential directions. We introduce the 
tangential gradient operator 6 = ,6n), where 6i = {5ij — ViUj)Dj. Applying this tangential 

operator to the boundary condition we have 

(2.4) DkudiVk + VkdiDkU = Spp, on dLi, 
hence we have 

(2.5) \Dt-i,u\ < C, on dLl, 
for any tangential vector field r. 

We next deduce the estimate for D^^^u on dLl. By a direct calculation, we have 

Lu = w^^{DijU — Dp^Aij{-,u,Du)Diu) 

= n- w'^^iAij - Dp^Aij{-,u,Du)Diu). 
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( 2 . 6 ) 




Differentiating the equation (|2.1I) with respect to we have, for A: = 1, • • • ,n, 

^ ^i^ij^k D]^Aij DzAijV,]^^ Dp^AijDiUji^ 

= + D^Buk + Dpi^BDiUk, 

which implies 

( 2 . 8 ) Luk = D^f^B + D^Buk + Dp^BDkiu + w'^^DkAij + D^AijUk), for k = !,■■■ ,n. 

If we consider the function h = UkDkU — ip{x,u), by (12. 6 p and (12. 8 p . we immediately have 

(2.9) \Lh\ < C(1 + ^ w** + iD^ul), in D. 

From the positivity of B we can estimate, 

(2.10) 1 < Cw^\ (wii)^ < Cw^\ 

Thus we obtain from (j2.9p and the boundary condition (|1.2n . 

( 2 . 11 ) \Lh\ < C{1 in Q, and h = 0 on dQ. 

From the uniform ^d-convexity of D (|1.12l) and the regularity of A , there exists a defining function, 
(f> G C‘^{Cl), satisfying (/) = 0 on dfl, Dcj) 7 ^ 0 on dQ. and (/> < 0 in D, together with the inequality 

( 2 . 12 ) Dijcj) - Dpi^Aij{-,u,Du)Dk(t) > 5il, 

in a neighbourhood M of dfl, whenever D^u > ip{x,u), where Ji is a positive constant and I denotes 
the identity matrix, with Af and di depending also on do, A and |u|i; 0 . We remark that ( 12 . 121 ) follows 
from (jl.l 2 l) . using the continuity of Dp A with respect to x and z together with an appropriate extension 
of the distance function, as in for example [also]. In particular, we can take (j) = —d + td? near 9D, for 
a large enough positive constant t, where d(x) = dist(x, dQ) is the distance function of D. Accordingly 

(2.13) L(j)>6iJ2w^\ 

for h > 0, d < do, for a positive constant do also depending on do) ^ and |u|i;n. By (12.lip . (I2.13P and 
choosing —(j) as a barrier function, a standard barrier argument leads to 

n —2 

D^h<C{l + M^-^), ondn, 

where M 2 = sup \ D'^u\, so that we have the estimate 
o 

72 — 2 

(2.14) < (7(1 + M 2 ) , on (9D. 

We conclude from (12.51) . (I2.14p and the ellipticity of u that 

72 — 2 

(2.15) < (7(1 + M 2 ) , on df}, 

for any direction We remark that if B is independent of p or n = 2, then the term in M 2 is not 
present in (I2.15p . 

We have now established the mixed tangential normal derivative bound and the double normal 
derivative bound on dQ so that it remains to bound the double tangential second derivatives on dQ. 
We shall adapt the delicate method in HZ], which is specific for the Neumann boundary value problem, 
to obtain the double tangential derivative bound on the boundary and consequently the global second 
derivative bound. 

Proof of Theorem \l.l[ First we note from the comparison principle. Lemma 3.1, that u > u in D or 

u — u is a constant. Discarding the second case, we modify the elliptic supersolution u by adding a 

perturbation function —acj) , where a is a small positive constant and (j) is the defining function of the 

domain D satisfying </) = 0 on dfd, (/> < 0 in D and Dyf) = —1 on 912. Note that the new function 
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u = u — acj) is still uniformly elliptic in 0 if a is sufficiently small. Also, by a direct computation, we 
have 

Dy{u — u) = DyU — DyU — aD^cj) 

(2.16) = - ip{-,u) + a 

> a, 

on where the non-decreasing of ip and u > u on dQ are used. If we define a function with the form 
<1> = with a positive constant K, we then have > Ka > 0 on dfl. We now introduce an 

auxiliary function v, given by 

(2.17) n = u(-,0 = - n'(-,0), 

for X G Q, 1^1 = 1, where a, k are positive constants to be determined, 

$ = 

ei 

is the barrier function in Lemma l2.II with the above constrncted u, and v' is defined by 


(2.18) = 2 (C • iy)C{Diip{-,u) - DkuDiUk - AijVj), 

with = ^ — (^ • Here is a C'^’^(il) extension of the inner unit normal vector field on dQ. The 
strategy of onr proof is to estimate u at a maximnm point in Cl and vector in the same form as 
(I2.15p . from this we conclude a corresponding global estimate for D^u in Cl from which follows the 
desired estimate (jl.lhp . 

Case 1. We snppose that v takes its maximum at an interior point xq G H and a vector C,- Let 

(2.19) 

then the function H also attains its maximum at the point xq € Cl and the unit vector C- The following 
analysis follows the method of Pogorelov type estimates in [30], with some modification, adapted from 
m, to handle the additional term v'. Accordingly we have, at the point xq. 


H = logx = log(rcgg — v') + ^\Du\‘^ + k<1>. 


Djiwcf — v') 

0 = DiH =-—— -h aDkuDikU + for i = 1 ■ ■ ■ n. 


( 2 . 20 ) 


0 > DijH = 


— v' 

- v') _ - v')Dj{w^^ - v') 

— v' — v'Y 

+a{DikuDjkU + DkuDijku) + nDij^, 


and consequently, at xq 
0>CH = 

( 2 . 21 ) 


1 


-L{w^^ - v') - Di{w^^ - v')Dj{w^^ - v') 


— v 

+aw^^ DikuDjkU + aDkuCuk + kL^. 

Next, we shall estimate each term on the right hand side of (|2.21ll . We start with some identities. By 
differentiation of the equation ( 12.111 in the direction we have in accordance with (| 2 . 8 p . 

w ^ DzAijU^ Dp^Aj^jDixi^') 

= D(^B + DzBu^ + Dp^BDiu^, 


( 2 . 22 ) 













(2.23) 


and a further differentiation in the direction of ^ yields, 

^[DijU^^ D^^Aij {^Dp^pj^Aij'jDf^u^Di'ii^ 

{D2;Aij)u^^ {^Dp^Aij'jDj^u^^ ^(^D^zAij'jii^ 

= w'-^w^'-D^WijD^Wki + + {DzzB){u^f + {Dp^p^B)DkU^Diu^ 

+2{D^^B)u^ + 2{D^p^B)DkU^ + 2(Dzpf,B){DkU^)u^ 

+{DzB)u^^ + {Dp^.B)DkU^^. 

Using (|2.23l) and the regular condition (11.31) . (see (3.9) in [3U]), we have 

(2.24) > w^^w^'-D^WijD^Wki - C[{1 + Wii)T + [waf], 

where we denote T = rc** to avoid any confusion with the usual summation convention. When calcu¬ 
lating LA^^, there will occur third derivative terms of u, which are controlled using (|2.22l) . We then 
obtain 

(2.25) \CA^^\ < C[{1 + w^^)T + Wii] 
and by a similar calculation, we have 

(2.26) \Cv'\ < C[(l + Wii)T + Wii], 

Combining (|2.24l) . (|2.25l) and (|2.26l) . we have 

(2.27) ^(w^CC - v') > w'-^w^^D^WijD^Wki - <^[(1 -h Wii)T -h (iCm)^]- 
By Cauchy’s inequality, we have 

(2.28) Di{w^^ — v')Dj{w^^ — v') < (1 -|- 9)w^^DiW^^DjW^^ + C{9)w^^Div'Djv' 
for any 9 > 0, where C{9) is a positive constant depending on 9. 

By (12.31) . (12.81) . (I2.27p and (I2.28p . we then obtain from (12.211) 

0 >- ^—^w^’^w^^D^WijD^Wki - 

1 


(2.29) 


— v‘ 

+awii + kT - C{ 
C{9) 


— V 
w^Wiv'Djv'. 


1 + I 

v'Y 

j[(l -h Wii)T + [wiif] +a + K} 


{w^i -v'Y 

Without loss of generality, we assume that {wij} is diagonal at xq with maximum eigenvalue ten. 
We can always assume that ron > 1 and is as large as we want; otherwise we are done. We proceed 
first to estimate the third derivative terms in (I2.29p . From the inequality (3.48) in [T7], we have 

(2.30) D^WijD^Wki - vA^DiW^^Djw^^ > 0. 

Wll 

Moreover since v' is bounded, ten and are comparable in the sense that for any 0 > 0, there exists 
a further constant C{9) such that 

(2.31) Itcii — -|-u'l < Stcii, 
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if wii > C{6). From (12.301) and (12.311) . we have 


(2.32) 


i-e 


D^WijD^Wki >- jw'^^DiW^^DjW^^. 


— V 


Next we use DiH = 0 in ()2.20p . to estimate 
(2.33) 


w^WiW^^DjW^^ < 2w^^[\Div'\'^ + (wgg - v')‘^{aDkuDikU + kA^)^] 
< C[w^^ + — v'Y{a^Wii + k^T)]. 


Using ()2.32l) and (I2.33p in (|2.29l) . together with (|2.31l) . we then obtain for wu > C{6), 
(2.34) awii + kT < C[a + k + (1 + a^9)wii + (1 + 


By choosing a, k large, and then fixing a small positive 9, we thus obtain an estimate Wii(xo) < C, 
which implies a corresponding estimate for |Zl^u(xo)|. 

Case 2. We consider the case xq € dQ, namely the function v(x,^) = — v') attains 

its maximum over fl at xq G 90 and a unit vector We then consider the following three subcases of 
different directions of For this we employ the key trick from m- 

Subcase (i). ^ = v, where v is normal to 90 at xq. Since from p2.14p we already obtained the double 
normal derivative bound, we have 

n —2 

(2.35) u(xo, z^) < (9(1 + M 2 ) , on 90. 


Subcase (ii). ^ is neither normal nor tangential to 90. The unit vector ^ can be written as 


(2.36) 


^ = (A r)T + (A 


where t € S"^ with r • = 0, (A t)'^ + (^ • ^)^ = 1 and ^ ■ u ^ 0. By the construction of v', we have 

clt Xq , 


(2.37) 


• T)^WrT + (? • + 2{C • r)(A v) 


Wr 


= (^ • rfwTT + (? • + v'{x, 0- 

By the constructions of u, we then have 


(2.38) 

which leads again to 


vixo,C) = (C • Tfv{xo, r) + (A t^fvixo, v) 
< {i ■ r)A(xo, 0 + (C ■ ^fv{xQ, v), 


v{xo,0 < u(xo, I') <C{1 + M 2 )'*-i , on 90. 
10 


(2.39) 



















Subcase (in). is tangential to d^l at xq. Observing the construction of v', we have v'{xo,^) = 0. 
We then have, at xq, 

0>D^v = - u')] 

= - v')Di,{^\Du\‘^ + k^>) + Dy{w^^ - u')] 

(2.40) = + ^0} 

= + aDku{(pk + <fzDkU - DiuDkVi)\w^^ + - Dy{A^^ + u')} 

> etl^“l^+'^'*’{(KCo - aM)w^^ + - Dy{A^^ + u')}, 

where cq = — , M = max \DkuUpk + (pzD^u — DiuDki'Al. The above inequality gives a relationship 
*^1 x£dn 

between w^^{xo) and D^u^^{xo), namely 

(2.41) < —(rcco — aM)w^^ + Dy{A^^ + v'), at xq. 

On the other hand, by tangentially differentiating the boundary condition twice, we obtain 


(2.42) Dku5i6jVk + 5iDku6jUk + 5jDku6iVk + VkSiSjDkU = 6i5jtf, on 512. 


Hence at xq, for the tangential direction ^ we have 

Di/U^^ ^ ‘l{5i'L'k)S)jku(^i(^j T C 

(2.43) > ip^Dijuiiij - 2{5iVk)Djkuiiij - C 

> ipzW^^ - 2{5iVk)Djkuiiij - C, at xq, 

where the double normal boundary estimate (I2.14p is used in the second inequality. The inequality 
(I2.43P clearly provides another relationship between DyU^^{xo) and w^^{xq). Combining this with 
(I2.4ip . we obtain 

(2.44) (kcq - aM + (pz)w^^ < 2{5iVk)Djku(,i(,j + D„{A^^ + v') + C, at xq. 


Without loss of generality, we can assume the normal at xq to be = (0, • • • , 0,1), and correspondingly 
we may assume {rcij(xo)}i,j<n is diagonal with maximum eigenvalue rcii(xo) > 1, as in the interior 
case. Observing that the first term on the right hand side of p2.44H only involves tangential second 
derivatives and using ()2.15p . we can then estimate at xq, 


(2.45) 


(kco — aM + (Pz)w^^ < C{wn + \DD^u\) 
< + eM2 + C*£. 

We now choose k sufficiently large, such that 


(2.46) 

and again we obtain 

(2.47) 


K > — [aM — inf (pz — C], 
Co 


c(xo,0 < <^(1 + M 2 )"- 


We now conclude from the above three subcases that if v attains its maximum over 12 at a point xq G 512, 

then v{xo,(,) is bounded from above as in (12.471) . which implies the second derivative D^^u{xo) is also 
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similarly bounded from above. Combining the above two cases, and using the Cauchy inequality, we 
obtain the desired estimate (jl.l5p and complete the proof of Theorem 11.11 □ 

As remarked in Section [U we can relax the supersolution hypothesis when Dp^A = 0, that is A is of 
the form (I1.17p . Moreover the details are then much simpler as we do not need to extend the Pogorelev 
argument to handle third derivatives. Here we proceed in accordance with Remark 1 in Section 3 of 
m, assuming as there initially that B is convex with respect to p, and replace the auxiliary function 
V in (j2.17p by 

(2.48) V = v(x,^) = — v' + — \Du\'^ + 

where now u € in <1> = is an elliptic function with ft > u in H, as in Lemma l2.II In 

place of ()2.27D . we now have the simpler inequality 

(2.49) L{w^^ - v') > -C{1 + T + Wii). 

And we obtain an estimate from above for if the maximum of v occurs at an interior point of H 
by taking again sufficiently large constants a and k. If the maximum of v occurs on the boundary 911, 
then we proceed as in Case 2 above except now the technical details are simpler and we do not need 
> 0 on 911 but we do need instead H uniformly convex or more generally (pz + 2ki > 0, where ki 
is the minimum curvature of 911, to use (12.4311 . We then obtain the estimate (11.1511 as before except 
that the dependence on u is replaced by a dependence on an elliptic function u. The removal of the 
condition that B is convex in p can then be addressed in the same way as in m by using Theorem 
O to construct a supersolution when B is replaced by its infimum and invoking the full strength of 
Theorem o 

Remark on Lemma 12.IL The proof of Lemma 12.11 following mm applies very generally. In fact, 
similarly to Theorem 2.1 in [23], we may replace the function “logdet” in (12.111 by any increasing 
concave function / on an open convex set T in the linear space of n x n symmetric matrices S”, 
which is closed under addition of the positive cone. Here the ellipticity conditions are replaced by the 
augmented Hessians Mu{Q)^ Mu[Cl) C T, which imply the operator F is elliptic with respect to u and 
ft on H and fl, respectively, and is replaced by T)-.. in the definition of L. The general case is 
covered with a slightly different proof in Section 4 of the forthcoming paper |8|; see also [TO] for the 
/c-Hessian case. However for the special case of m, the proof of Lemma l2. II from [7119| may also be 
simplified somewhat by avoiding the perturbation of u that is one of the key ingredients of the general 
argument used there. To see this, we may modify the calculations in the proof of Lemma 2.2 in [7], 
with e = 0 and v = u — u, (without using concavity!), to arrive at the inequality, 

(2.50) > Ke^^{w^^[DijU - Aij{-,u,Du) - Wij] - - Dp^B{■,u, Du)Div}, 

for any positive constant rj and sufficiently large constant K depending also on rj. We then obtain (12.3p 
using the simple inequality 

w'^^[DijU — Aij{-,u,Du)] > w^^X[Mu] > 0, 

where X[Mu] denotes the minimum eigenvalue of Mu, and taking p sufficiently small. 

3. Existence and solution estimates 

In this section we complete the proof of Theorem 11.21 and provide alternative conditions for the 
maximum modulus for solutions of the Neumann problem (ll.ip - (ll.2l) . First we formulate a comparison 
principle for general oblique boundary value problems (|1.7p - (|1.8p with F defined by (11.91) . with A and 
B non-decreasing in z, and G € (7^(911 x M x M”, non-increasing in 2 . 
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Lemma 3.1. Let u,v € C‘^{Ll) PI (7^(0) with T elliptic, with respect to u, in Q and Q oblique with 
respect to [u, v] on dQ, where [u, ?;] = {9u + (1 — 0)v : 0 < 0 < 1}. Assume also that either G is strictly 
decreasing in z or A or B are strictly increasing in z. Then if iF[u] > T-[v\ on the subset of LI where 
F is elliptic with respect to v and Q[u\ > Q[v\ on dLl, we have 

(3.1) u < V, in ri. 

Moreover if we assume that F is elliptic with respect to [u, v] on all of LI , we may relax the strict 
monotonicity condition on A, B or G, provided u — v is not a constant. 

The proof of Lemma 13.11 is standard. By approximating by a subdomain and approximating u 
by a smaller elliptic function u satisfying F\u\ > F\v\, we infer that the function u — v can only take 
a positive maximum on the boundary dLl and m then follows from the obliqueness and the strict 
monotonicity of G. When G is only non-increasing in z, then we can take u = u — e{(f — minc/i) for a 
dehning function cf G G‘^{Ll) nC'^(fi) such that (f = 0 on dLl, (/> < 0 in and sufficiently small e > 0, to 
ensure G[u\ > G[u] on dLl, whence a positive maximum of u — v must be taken on in LI and we conclude 
mi from the strict monotonicity of F with respect to 2 ;. Note that when G is strictly decreasing, we 
need only assume G is weakly oblique, that is Gp • z/ > 0 on dLl while when F is strictly decreasing 
we need only assume F is degenerate elliptic. In the case when there is no strict monotonicity, the 
difference w = u — v will satisfy a linear uniformly elliptic differential inequality of the form 

Cw := DijW + biDiW + cw >0, 

together with an oblique boundary inequality, /3 • Dw > 'jw, with coefficients c < 0 and 7 > 0, and the 
result follows from the strong maximum principle and Hopf boundary point lemma; (see [5]). 

From Lemma [ 3 .II we have immediately the uniqueness in Theorem ll.2l and the inequality u<u<u, 
where u and u are the assumed elliptic supersolution (I1.13l) - (|1.14p and subsolution. 

Next we obtain a gradient bound for A-convex functions for Neumann problem (ll.ip - (ll.2p . where 
A satisfies a quadratic bound from below, (|1.18p . by a modification of our argument for the Dirichlet 
problem in [S]. For this purpose, we formulate the following gradient estimate as a lemma. 

Lemma 3.2. Let u G G‘^{Ll) n C^(f2) satisfy 

(3.2) D‘^u > -/io(l + \Du\^)I, 
in a G‘^ domain LI C M"’, with 

(3.3) > —a, 

on dLl, where /tq and a are non-negative constants. Then we have the estimate 

(3.4) \Du\ < G, 
where C depends on fVo,a,Ll and sup|n|. 

Proof. Defining u = u — acf, where as in Section [21 (/) G G^{Ll) is a negative dehning function for Ll 
satisfying = — 1 on dLl, we see that v ■ Du > 0 on dLl. Consequently at a maximum point xq G D 
of the function 

(3.5) w = e^^\Du\^, 
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we have 


(3.6) Du-Dw< 0. 

From (13.2p . we have 

D'^u = D‘^u — aD'^cj) 

.3^^ >-no{l + \Du\‘^)I -aA^I 

>-Ho{l + 2\Du\‘^ + 2a^\D(l)\‘^)I-aA^I 
> + |F>u|^)/, 

for some positive constant fj,i depending on fiQ, a, Dcj) and A,^, where A^p denotes the maximum 
eigenvalue of the Hessian matrix of (j) and depends on the domain H. With the lower quadratic bound 
(13.71) for the Hessian matrix D^u in hand, by choosing the constant k sufficiently large as in Section 4, 
[S], we can obtain from (13.61) . 

(3.8) \Du\ < C, 


at xo, where the constant C depends on ^q, a and H. We then conclude a global gradient estimate 
from (13.81) and the construction of u, 

(3.9) |H>m| < C, 

where C depends on /io, H, a and sup |m|. □ 

We remark that by taking more careful account of the constant dependence in the proof of Lemma 
13.21 we infer a sharper estimate 

(3.10) \Du\<C{l + cr), 
where C depends on /tq, H and sup |n|. 

Note that the gradient estimate (|3.4p in Lemma 13.21 and the sharper gradient estimate (j3.10l) hold 
for any solution u satisfying the weak convexity condition (|3.2I) and the lower bound condition (13.31) 
for normal derivative on the boundary. We now apply Lemma 13.21 to obtain the gradient estimate for 
A-convex solutions of the Neumann problem (in])-([L2l) with A satisfying the lower quadratic bound 
(jl.181) . From the A-convexity of the solution u and the quadratic structure condition (jl.l8p , the solution 
u satisfies the weak convexity condition (j3.2p . The Neumann boundary condition (IE21) provides us a 
lower bound D^u > ip{x,u). Applying Lemma [321 we then obtain the global gradient estimate 
for Neumann problem (|l.ip - (|1.2l) . that is \Du\ < C for C depending on /tq, H, cp and sup |ri|. 

Since we now have obtained the derivative estimates up to second order, we can use the continuity 
method to prove our existence theorem. 


Proof of Theorem \1.2i . From the second derivative estimate, Theorem 11.11 and the preceding solution 
and gradient estimates we can derive a global second derivative Holder estimate 

(3.11) |n|2,«;n < C, 

for elliptic solutions u G (7^(17) nC^(H) of the semilinear Neumann boundary value problem (ll.ll) - (ll.2ll 
for 0 < a < 1. The estimate (13.lip is obtained in [16], Theorem 3.2, (see also [T5l|23|). With this 
estimate, we can use the method of continuity, Theorem 17.22 and Theorem 17.28 in [5], to derive the 
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existence of a solution u € using the supersolution u as an initial solution. To be rigorous, we 

should assume that A and B are (7^’“ smooth, (p is smooth and Vl £ for some a > 0 to get 
a solution u G by the Schauder theory, (see [5], Section 6.7), and then by approximation get a 

solution u £ Alternatively we can use the Aleksandrov-Bakel’man maximum principles (see 

[5], Theorem 9.1, Theorem 9.6) to carry over the proof of Theorem 1 1.1 1 to solutions u £ lT^’”(n)nC^(ii) 
and use Lp regularity as well, Section 9.5) to improve C'^’"(n) solutions with 0 < a < 1 to be in 
the Sobolev spaces D C'^’^(fi) for all p < oo, 0 < J < 1. □ 

In the rest of this section we will consider more explicit conditions for solution bounds. Here we 
consider the oblique boundary value problems ()1.7l) - (|1.8p with F defined by (|1.9p and G defined by 
(I1.16p . that is the Monge-Ampere type equation (jl.ip together with the oblique boundary condition 

(3.12) Dpu = ip{x,u), on dQ. 

First we note that we also obtain bounds for solutions u of (ini-disi) if u and u are only assumed to 
be supersolutions and subsolutions, without any assumed boundary conditions, provided we strengthen 
the monotonicity of p. In particular we may assume, as in m, there exists a positive constant 70 such 
that 

(3.13) Pz(x,z)>'yo 

for all (x, z) £ dQ x M. In the light of Lemma l3.ll we may interpret a supersolution as satisfying (11.131) 
only at points of ellipticity. Since A and B are non-decreasing, supersolutions and elliptic subsolu¬ 
tions are preserved under addition and subtraction respectively of positive constants. Accordingly, by 
subtracting a positive constant from u and using (|3.13p we can assume Dpu > p{x,u) on dCl, whence 
u > u in H. Similarly by adding a positive constant to u we obtain Dpu < ip{x, u) on dQ, so that u <u 
in H. Note that for this argument we may replace (13.131) by the weaker conditions 

(3.14) (sign2:)(^(-, z) —>■ 00, as \z\ —>■ 00. 

The conditions (|3.13p . (|3.14l) may be further weakened when constants are subsolutions or supersolu¬ 
tions. We first consider the bound from below, under the following conditions: 

(3.15) A(x,z,0) < 0, det[—A(x, z, 0)] > B{x,z,0), for all x £ 11, z < —K, 

(3.16) ip{x,z) < 0, for all x £ 911, z < —K, 

where 77 is a positive constant. Under the assumptions (13.151) and (13.161) . we can readily obtain the 
solution bound as follows. Suppose u attains its minimum over H at a point xq and u{xq) < —K. If 
xo £ U, we have Du{xq) = 0, D‘^u{xo) > 0. From the equation (II.ip . we have det[—A(xo, m(xo), 0)] — 
B{xo,u{xo),0) < 0 so that by (13.151) . we must have m(xo) > —K. If xq £ 911, we have Dpu{xo) > 0. 
From the oblique boundary condition, we have p{xo, m(xo)) > 0. By (I3.16p . we again have u{xq) > —K. 
Note that condition (j3.15l) implies sufficiently small constants are subsolutions of the oblique bound¬ 
ary value problem (jl.ll) - (|3.12p thereby providing lower solution bounds, by the comparison principle. 
Lemma 13.11 Therefore the subsolution assumption in Theorem 11.21 can be replaced by the structure 
conditions (j3.15p and (I3.16p . with minu replaced by —K in X. We also remark that condition (|3.15p 
follows from a uniform monotonicity condition on A, namely 

(3.17) D^Aij{x,z,p)(iCj >7i|^P, 

for all {x,z,p) £ U X M X R”, ^ £ R'^ and some 71 > 0, which is a stronger form of the A4w condition 
used for generated prescribed Jacobian equations in geometric optics in [3 EH], together with B being 
non-decreasing in z. 

In this sense, the condition (13.151) is a weakening of the uniform monotonicity of A, while the 
condition (|3.16p is a weakening of the uniform monotonicity of p. On the other hand, condition 
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(j3.15p is restrictive in that it excludes the case when A is independent of z, which occurs in optimal 
transportation. 

Corresponding conditions also provide bounds from above. Here though the analogue of (|3.15p is 
more general, namely 

(3.18) det[—H(x, z, 0)] < B(x,z,0), for all x € H, z > K, A(x,z,0) < 0, 
while instead of (j3.16D . we have 

(3.19) (p(x, z) > 0, for all x E dO., z > K, 

where K is a positive constant. Note that condition (|3.18p extends the condition in Section 4 of [9], 
namely that the maximum eigenvalue of A(x, z, 0) is non-negative for all x & Q, z > K for some positive 
constant K and implies that constants larger than K will be supersolutions, where they are elliptic. 

To complete this section, we derive a lower bound for optimal transportation equations and present 
the corresponding existence result. 


Optimal transportation equations. In the optimal transportation case, we can replace the existence 
of a subsolution in Theorem 11.21 bv an extension of the sharp conditions (1.4), (1.5) in [T7], through 
an extension of the Aleksandrov-Bakel’man estimate in Theorem 2.1 of m- Optimal transportation 
equations are special cases of prescribed Jacobian equations where the mapping Y is generated by a 
cost function c defined on a domain V C M” x M"'. We assume Cl x A G V, for some domain A C M”, 
and c E C‘^{V) satisfies the conditions, (from [22]): 

Al: For each x E H, the mapping Cx{x, •) is one-to-one in y E P* = {y E ]R”'|(x,y) E P}; 

A2: detCa;,y / 0 on P. 

Then the mapping Y is given by 

(3.20) Y{x,p) = c~^{x,-){p) 

and is well defined for p E Pa, = {p G 1K”'| P = Cx{x, y) for some y E P^}. In the resultant Monge-Ampere 
type equation, we then have from (ll.20|l . 

(3.21) A{x,z,p) = A{x,p) = Cxx{x.,Y{x,p)), B = \(leicx,y\'4), 

and equation dni is well defined for solutions u which are A-convex and satisfy Du{x) E lAx, for each 
X E H. We call such solutions admissible. In the optimal transportation case, c-afhne functions, that is 
functions of the form u = c(x, y) J-cq, for constant cq and (H, {y}) C P are automatically supersolutions 
as they satisfy the homogeneous equation 

(3.22) det(P^u — A(x, Pu)) = 0, 


and hence provide upper bounds for solutions of (weakly) oblique boundary value problems, 

(3.23) D/su = (p{x,u), on dQ, 

where ■ v >0 on dil, under a uniform monotonicity condition (I3.13p . For lower bounds we impose a 
structure condition 

(3.24) 


ipix,z,p) < 


fix) 


f*oY{x,p) 

for all X E H, z < mo, Y{x,p) E A, where / > 0, E L^(H), /* > 0, E L\^^{A) satisfy 


(3.25) 


.[^<1 


r 


and mo is a constant. 

We now have the lower solution bound in the optimal transportation case. 
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Lemma 3.3. Let u € C‘^{Ll) H be an admissible solution of equation (I1.19p . in the optimal 

transportation case (I3.20p . with cost function c satisfying Al, A2. Suppose that ip satisfies (|3.24p and 


(3.26) Dgu < 7on + ipo on dLl, 

for u < niQ, where (5 G L°°{dLl), jd ■ n > 0 on dLl and 70 > 0 and (po >0 are constants. Then we have 
the lower bound 

(3.27) u > -C, in n, 

where C is a positive constant depending on 0,/,/*,/3,70 , (/Jq o,nd c. 


Proof. Our proof is adapted from the second author’s 2004 Singapore Institute of Mathematical Sciences 
lectures and the case where c(x, y) = x ■ y, that is Y = p and ^ = 0, in |17] . First, we note that if we 
have a global support from below at a point xq G Ll, that is 


(3.28) 


u{x) > u{xo) + c(x, yo) - c(xo,yo) 


for all X G n, then we must have yo = Y (xq, Du{xo)). Dehning T = Y{-, Du), we have by (|1.19D . (13.2411 
and the change of variable formula 


(3.29) 



det DT\f*oT 


> / r 


JT{no) 

where Hq = {a; £ 'u{x) < mo}. Hence by our condition (|3.25l) on / and /*, there exists a 

yo G A — T{Llo). It then follows by upward vertical translation of a c-affine lower bound, that 
exists a point xq G dLlo such that 


point 

there 


(3.30) u{x) > u{xo) + c(x, yo) - c(xo, yo) 

for all X G II. If xo G dLl, we must also have 


(3.31) 


Dpu{xo) > Di3c{xo,yo) 


whence by the boundary inequality (|3. 2611 .we obtain 


(3.32) 


u{xo) > —[Diic{xo,yo) - To]- 
70 


If Xo 0 dfl, then we must have u(xo) = mo. Hence by (13.3011 again, we obtain for xq G dQ 


(3.33) 


u(x) > u(xo) + c(x, yo) - c(xo, yo) 

> — lDfjc(xo, yo) - < 7 * 0 ] + c(x, yo) - c(xo, yo) 
70 

> - (— +diamII)sup|c3;(-,yo)| 

7o 7o n 


while for xo 0 <911 we obtain 
(3.34) 


u(x) > mo — diamllsup |ca;(-, yo)|. 

fi 
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To remove the dependence on yo in (|3.33l) and (I3.34p . we may consider an exhaustion of A, say by 
defining subdomains 

(3.35) Ar = {y € A\ \y\ < R,dist{y,dA) > 
for i? > 1. Then by (I3.25p . we have 

(3.36) [ f= [ f* 

for some sufficiently large -R, and we obtain from (I3.33D and (|3.34|) . the estimate, 

(3.37) u{x) > min{mo, — —} — (— + diamhl) sup |i7c|. 

7o 7o nxAfl 

This completes the proof of Lemma 13.31 □ 

As a corollary of Lemma 13.31 and the proof of Theorem 11.21 we then have the following variant of 
Theorem ll.2l in the optimal transportation case. For this purpose we note that the boundary condition 
(jL2p and the monotonicity condition (j3.13p imply ()3.26l) with /3 = and 


7^0 = -7oW'0 + sup (/?(•, mo). 

an 

Corollary 3.1. Suppose that equation (11.111 is a prescribed Jacobian equation of the form (11.1911 gen¬ 
erated by a cost function c € satisfying conditions Al and A2 and lAx = K"" for all x G 0, with 

satisfying the structure conditions (I3.24p . (|3.25p . Let A, B, ip and 12 satisfy the hypotheses of Theorem 
11.ill except for the existence of an elliptic subsolution, with p satisfying (13.1311 and 12 assumed to be 
uniformly A-convex with respect to p and —C, that is (11.1111 holds for p ■ u > p(-, —C) on dQ, where 
C is the constant in Lemma 3.3. Then the Neumann boundary ualue problem dnp-don has a unique 
elliptic solution u G C^’"(fi) for any a <1. 

We remark that as in [m, condition (13.2511 is necessary for an elliptic solution u G (7^(12) n C‘^’^(12) 
of (ll.iyp . with Du{x) G Ux for all x G 12. 

In accordance with our remarks following the statement of Theorem 11.21 pertaining to the special 
case (jl.l7ll . and using the argument at the end of Section [2l we can remove the supersolution condition 
in Corollary 13.11 for convex domains. To apply the argument at the end of Section [21 we also need to 
use the existence of an elliptic function, as provided by Lemma 2.1 in [7]. In this way, we obtain an 
extension of Theorem 1.1 in which corresponds to the special case c(x, y) = x - y, (or equivalently, 
the case c{x,y) = —|x — yp/2). Note that the matrix A generated by the cost function satisfies (I1.17P 
when the cost c = c(x — y). Examples of regular and strictly regular cost functions are given in |30j 
and m- However most of these examples do not satisfy ILx = M” and in general we need additional 
controls on gradients to prove classical existence theorems. 

We also remark that Lemma 13.31 and Corollary 13.11 are readily extended to generated prescribed 
Jacobian equations [28] . 


4. Oblique boundary value problems 


In this section we consider more general oblique boundary value problems for Monge-Ampere type 
equations under the hypothesis that the matrix function A is strictly regular. As remarked in Section 
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dl this condition also leads to a much simpler proof in the Neumann case. Also we do not need to 
restrict to semilinear problems of the form (jl.l6|l but can consider nonlinear boundary conditions of 
the general form (|1.8h . where G is also concave with respect to p. Our approach is already indicated 
in Section 4 of [30] and we will carry over some of the basic details from there. Moreover our results 
can also be seen as special cases of those for general augmented Hessian equations in |8]. For second 
derivative estimates, we will assume that the function G G G^{d^l x M x M”) is oblique with respect to 
a solution u, that is from (|1.22l) . 

(4.1) Gp{-,u, Du) ■ v > Pq, on dil, 

for a positive constant /3o, and is concave in p, with respect to u, in the sense that 

(4.2) Gpp{-,u, Du) < 0, on dQ. 

We now have the following extension and improvement of Theorem 11.11 in the strictly regular case. 

Theorem 4.1. Let u G G^{Q) PI (7^(0) he an elliptic solution of the boundary value problem a- 
in a G^’^ domain Q C M”, which is uniformly A-convex with respect to G and u, where A G 
C^(H X M X M"") is strictly regular in Cl, H > 0, G (7^(0 x M x M") and G G G^'^{dLl x M x R”) satisfies 
mi) and dMI). Then we have the estimate 

(4.3) sup \D‘^u\ < G, 
where G is a constant depending on n,A,B,G,Ll,l3Q and |u|i;Q. 


Proof. As in the proof of Theorem 11.11 we first consider the estimation of the nontangential second 
derivatives. In the semilinear case (jl.lbp . we can simply replace n hy jS there and deduce in place of 
(I2.15D . the estimate 

71 — 2 

(4.4) < (7(1 + M2) "-1, on (9H, 

for any direction where as in Section O M 2 = sup^^ \D'^u\. In the general case, we have the same 
estimate m, from the estimate (4.4) in [30] . where now fd = Gp{-,u, Du). Now differentiating the 
boundary condition (USD twice with respect to a tangential G^ vector field r we obtain as in the 
estimate (4.10) in [HU] . 

/. '^TTjS P Dpi_p^GukT-uij- f7(l + M 2 ) 

>—(7(1 +M 2 ), on cin, 

by virtue of the concavity of G with respect to p. For convenience we write here Uir = UijTj, Utt = 
UijTiTj, Urrg = UijkTiTjfdk- To handle the pure tangential derivatives we extend the vector held r 
to all of 11 and set 


(4.6) V = Wrr - K{1 +M2)4>, 

where as in the proof of Theorem 11.11 cp G G‘^{Cl) is a negative dehning function for 11 satisfying 
Di^(f) = — 1 on dLl and AT is a constant such that 


(4.7) 


DjsiwijTiTj] > -K{1 + M2)/3 o, on dLl. 
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In particular we may fix r with Ti = Xi — {x ■ v)vi, i = 1,..., n, where as in Section [2l i/ is a smooth 
extension of the inner normal v to 0. It then follows that D^v > 0 on dVL so that v must take its 
maximum on 11 at an interior point xq € with Cv{xq) < 0. Now we can adapt the proof of the 
interior second derivative estimate in |22j and [29], differentiating the equation dni), in the form (I2.1li . 
twice with respect to r and using also the concavity of the function “logdet”, together with (14.51) to 
control K, to estimate at xq, 

(4.8) Aij^kiUkrUiT < <^[(1 + + |DurP]- 

We note that when we twice differentiate (|l.ip with respect to a variable vector field r, to calculate 
Cv, we encounter terms arising from derivatives of r which are not present in the constant case (12.231) . 
Apart from the terms in third derivatives these can be directly estimated by (7(1 + M 2 )rc**. Retaining 
the third derivative terms, we would supplement the right hand side of (14.81) . by 

DrWijDrWkl + Aw’^WiTkDr'Wjk 
= -w'^^wADrWijDrWkl + ^w'^^W^'-WjkDiTkDrWkl 

< -w'^^W^\DrWij - 2WjkDiTk){DrWkl - 2WilDjTl) + Aw^’^W^'-WilDjTlWjkDiTk 

< A{DiTif 

so that the estimate (14.81) is unaffected. To use the strictly regular condition, 

(4.10) Aij^kiiiijVm > 


for all ^,rj £ satisfying the orthogonality ^ T ry, where cq is a positive constant depending on A and 
we choose coordinates so that w is diagonalised at xq, so that 

w'-^ Aij^klWkrWlr = W^^Aii^klWkkWllTkTi 

(4.11) SE w'^'^Aii^klWkkWllTkTl - CM 2 

k,Di 

> cqu;** '^{wkkTkf - CM 2 . 

Hence we obtain from (|4.8I) . (14.111) and (|2.11l) . 

(4.12) DrrU{xo) < (7(1 + M2)2. 


At this point we need to return to our choice of (j) to ensure that inf 4> > —e for some small positive 
constant e. This can be done for example by mollification of the function — inf{d, e} for sufficiently 
small e, where the constant (7 = (7^ in (|4.12p will depend also on e. Alternatively, we may simply 
restrict to a boundary strip = {(/> > —e} and use the interior second derivative estimates |22l |30] to 
estimate v on the inner boundary {0 = —e}. Accordingly we obtain from (|4.12l) . 

v{xo) < Ce{l + M 2)2 + eM 2 


DrrU < (7e(l + M 2)2 + eM 2 on dQ. 
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and hence we get an estimate 
(4.13) 












Since for any direction we have, 

(4.14) = Urr + b{Url3 + UfSr) + 
where 

b=k-^, r = C-bl3, 

P ■ 7 

we then obtain a boundary estimate in the form, 

(4.15) sup |Zl^u| < eM 2 + C^, 

dn 

for any sufficiently small e > 0, by combining (|4.4I) and (14.1311 . The global second derivative estimate 
now follows from the global second derivative estimates in |29[ [30] by choosing e sufficiently 
small. □ 

The details in the proof of Theorem 14.11 can be further varied. For example we can replace v by 

(4.16) V = {1 — K(j))WrT, 

for a sufficiently large constant iF, where 4> is the same negative defining function as in the proof of 
Theorem El As remarked in Section [H we also obtain a much simpler proof of Theorem 11.11 in the 
strictly regular case, without need for the supersolution and monotonicity hypotheses. Moreover by 
flattening the boundary dfl in a neighbourhood A/" of a fixed point xi € cAl, we can localise the second 
derivative estimate by modifying (|4.16l) 

(4.17) V = r]{l — K(j))wrT 

where r/ is a suitable cut-off function satisfying Dj^t] = 0 on A/" (7 dQ. Accordingly, we obtain for any 
ball B = Bji{xo) of radius R > 0 and centre xq, the local estimate 

(4.18) \D^u{xo)\ < 

for elliptic solutions u € C'^{B D 11) n C^{B n Q) of (|l.lh satisfying (|1.2p on (7 <911, where B 7 511 is 
uniformly A-convex with respect to G and u in the sense that 

{DiUj - Dp^Aij{-,u, Du)uk)riTj < -Jq 

on i? 7 dQ for G{x,u, Du) > 0 and any unit tangential vector r and a positive constant 6o- The 
constant C in (I4.18P depends on n, A, B,Q,,So, cp and |u|i;Q. We also point out that comparability of 
differentiation with respect to a general vector field and a constant vector field in the proof of 14.11 
which follows from the identity (14.91) . is special to the Monge-Ampere case. A different and more 
detailed proof of the critical tangential estimate (|4.13l) is provided for more general augmented Hessian 
equations in [8|, Lemma 2.3. 

Returning to the example from conformal geometry in Section (T] namely (11.2411 , (11.2511 with A4 = 
H C M”, the A-convexity condition also simplifies in that H is uniformly A-convex with respect to G 
and u if and only if 

(4.19) Ki >—ce~'^ + hdn on 

where ni denotes the minimum curvature of clH, and Theorem 14.11 extends the second derivative 
estimates in m for this special case with c > 0. We remark though that the strictly regular case in 
Theorem 14. II also extends to general augmented Hessian equations and corresponding second derivative 
estimates for ()1.24ll for general / are proved in [8]. 

From Theorem 14.11 we can obtain existence theorems, which also extend Theorem 11.21 and Corollary 
13.11 in the strictly regular case. First we prove an appropriate extension of the gradient bound Lemma 

E 
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Lemma 4.1. Let u G C‘^{Lt) satisfy (13.2p in a domain Q, C M" and 

(4.20) \D/3u\ < (To, /3 • z/ > /3o 

on dQ, where /3 € L°°{dLl), |/3| = 1 and do and /3o are positive constants. Then we have the estimate 

(4.21) \Du\ < C, 

where C depends on and sup|n|. 

Proof. Invoking the tangential gradient 6u, we have the formula 

(4.22) DiyU = — — {Dpu — P ■ 6u) 
so that we can estimate 

(4.23) |L)n| < ^(|(5n| + do) + |(5n| 

Po 

on dLl, whence from (13.2p . we obtain 

(4.24) D‘^u> + 

on dQ, for a further constant pi, depending on po, /3o and do. Now we consider in place of (13.5p . the 
function 

(4.25) w = e™|(5np, 
so that at a point xo G dLl where w is maximised we have 


0 = 6u ■ 6w 

= (k\6u\‘^ + 2diu6juSiSju) 

4.26 ^ ' ^0^3) 

= e^'^[K\6u\^ + 25iu5ju{DijU — DyuSiVj)] 

> e'^'^{K\5u\^ — 2pi\5u\^{l + |(5np) — Cldnp], 

from P4.22I) and (I4.24p . where C is a constant depending on /3o, do and dLl. By choosing k sufficiently 
large we conclude the estimate p4.2ip on dLl and the estimate in all of Ll then follows from or Lemma 

[321 □ 


Lemma l4 .1 1 provides an extension of Theorem 2.2 in [T7| to the weaker convexity condition p3.2p . If 
we assume a stronger quadratic control from below on the Hessian, namely 

(4.27) Dijuiiij > -po{l + [L'^ap) 

for some constant pQ and any unit vector we can reduce to Theorem 2.2 and the corresponding 
remark in m as condition (14.271) implies that the function is semi-convex for large k. We also 
remark that the gradient estimates in Lemma 13.21 and Lemma l4.II have local versions. In particular, if 
we fix any ball B = Bji{xo) of radius R and centre xq G Cl, and suppose u G C'^{Q D H) D (7^(11 D B) 
satisfies (|3.2I) in H D H and (14.201) in dLl n B, then we have an estimate 

\Du{xo)\ < 


(4.28) 


22 













where C depends on fiQ,ao, (3o,Q and sup |u|. To prove (14.281) we modify our proof of the global estimate 
Lemma l4.1l bv maximizing in place of the auxiliary functions in [S] and (|4.25l) above, the functions 

(4.29) wi = W2 = 

over n n i?, dQ H B respectively, where rj G Cq{B) is a cut-off function chosen so that 0 < r/ < 1, 
r/(xo) = 1 and \Dri\ < 2/R. 

Note that (14.271) is satisfied in the special case (ll.25p so we obtain, for solutions of (11.241) . (11.251) . 
both local and global, gradient and second derivative estimates in terms of fl, and sup |u|. 

In order to apply Lemma l4. 11 we also need to assume that G is uniformly oblique in the sense that 

(4.30) Gp{x, z,p) ■ ly > l3o, \Gp{x, z,p)\ < ao on dCl, 

for all X € Q, \z\ < Mq, p € M"’ and positive constants /3o and do, depending on the constant Mq. Using 
the mean value theorem, we can thus write G in the semilinear form (|1.16l) so that Lemma l4. 11 as well 
as the solution estimates in Section [3l are applicable. 

We then have the following analogue of Theorem 11.21 with a much weaker supersolution condition. 

Theorem 4.2. Suppose that A, B, G and 11 satisfy the hypotheses of Theorem \4-l\ with G uniformly 
oblique satisfying (|4.30l) and concave in p for all (x, z,p) € dTl x R x R"’. Assume also that A and B are 
non-decreasing in z, G is strictly decreasing in z, A satisfies (ll.lSp and that there exists a supersolution 
u and an elliptic subsolution u of equation dEU) in G‘^{Tl) n C^(n) satisfying Q[u] < 0 and > 0 
respectively on dO. with H uniformly A-convex with respect to G and X = [u,u\. Then the boundary 
value problem dnp-dis]) has a unique elliptic solution u € C'^’“(n) for any a < 1. 

Analogously to Corollary 13.11 we also have from Lemma 13.31 an existence theorem in the optimal 
transportation case. Here we may also extend the condition (|3.13p by assuming there exists a positive 
constant 70 such that 

(4.31) Gz{x,z,p) <-Jo 

for all (x, z,p) G cAl x R x R”'. 

Corollary 4.1. Suppose that equation (11.11) is a prescribed Jacobian equation of the form (|1.19p gen¬ 
erated by a cost function c G G‘^{V) satisfying conditions Al and A2 and Ux = R*^ for all x G 11, 
with if satisfying the structure conditions (13.241) . ()3.25p . Suppose also that A,B,G and H satisfy the 
hypotheses of Theorem \4-1\ with G uniformly oblique satisfying (|4.30p . uniformly monotone satisfying 
(I4.3ip and concave in p for all {x,z,p) G 911 x R x R”, A satisfying (jl.lSp . B non-decreasing and 
H uniformly A-convex with respect to G and —G, where C is the constant in Lemma \S.A Then the 
boundary value problem (fLTD-lfLSl) has a unique elliptic solution u G C'^’“(ll) for any a < 1. 

Finally we remark that when G is assumed uniformly concave with respect to p, we only need A 
to be regular in Theorems 14.1114.21 and Corollary 14.11 and the global second derivative estimates follow 
exactly as in Section 4 of [3D]; see also [33]. Also the proof of Theorem 14.11 would carry over to the cases 
when G is non-increasing and A is non-decreasing, with either G^ sufficiently small or DzA sufficiently 
large and A again only assumed regular, (using in the first case the existence of an elliptic function 
and Lemma l2.ip . 
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